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Abstract. We present ISO observations of the S(0), S(l), 
S(2), and S(3) rotational lines of molecular hydrogen from 
two molecular clouds near the Galactic Center (GC). We 
have also measured continuum dust emission at infrared 
wavelengths with ISO and the rotational radio lines J=l-0 
of 13 C0 and C 18 and J=2-l of C 18 with the IRAM- 
30m telescope. Using the dust continuum spectra and the 
CO lines we derive a total visual extinction of ~ 15-20 
magnitudes toward these GC clouds. After correcting the 
H2 data for extinction, the gas temperatures are ~ 250 K 
and the column densities of warm gas are ~ 2 x 10 21 cm~ 2 . 
This is the first direct measure of the H 2 column densities 
of the warm component; with this, we estimate an NH 3 
abundance in the warm gas of - 2 1(T 7 . The column 
density of warm gas is, at least, a factor of 100 larger 
than the corresponding column densities derived from the 
warm dust. The observed ortho-to-para ratio (OTPR) is 
~ 1, clearly below the local thermodynamical equilibrium 
(LTE) OTPR for gas at 250 K of - 3. Low velocity shocks 
(~ 10 km s _1 ) are the most likely explanation for the 
column densities of warm gas and dust and the non-LTE 
H 2 OTPR. 

Key words: ISM: clouds - ISM: molecules - ISM: dust, 
extinction - Galaxy: center - Infrared: ISM: continuum - 
Infrared: ISM: lines and bands 



1. Introduction 

The central ~ 6° of our galaxy exhibit a large accumula- 
tion of molecular material which is forming big molecular 



clouds whose masses and sizes are so large as 10 6 M Q and 
15 pc, respectively. These clouds are denser (average den- 
sities of 10 4 cm~ 3 ), more turbulent (line widths of ~ 20 
km s _1 ), and hotter (with a warm component with tem- 
peratures, T, up to 200-300 K) than the c louds of the disk 
of the galaxy (see e.g. Morris & Serabyn 1996). The high 
temperatures in Galactic Center (GC) clouds were known 
basically by observations of NH3 inversion lines over lim- 
ited regions (Glisten et al. [1981 ; Mauersberger et al. 1986 ). 
Hiittemeister et al. (1993) analyzed 36 molecular clouds 
distributed all along the Central Molecular Zone and the 
"Clump 2" complex; they showed that high kinetic tem- 
peratures are a general characteristic of the GC clouds and 
not only of those located close to Sgr A and Sgr B2. In the 
disk of the galaxy, kinetic temperatures higher than 100 
K are associated with infrared sources, that is, embedded 
stars which heat the dust and subsequently the gas by 
collisions with the dust grains. The typical sizes of such 
regions are less than 1 pc. The high kinetic temperatures 
in the GC clouds are found in regions of ~ 10 pc, where 
one measures large col umn d ensities of cold d ust (T < 30 
K, Odenwald & Fazio |l984[ Cox & Laureijs |198S| ). This 
rules out gas-dust collisions as a possible heating mech- 
anism of the warm component. Dissipation of turbulence 
due to shocks induced by the rotation of the galaxy could 
be the main heating mechanism in the GC clouds (Wilson 
et al. |1982|). 



Unfortunately the NH3 abundance in the warm com- 
ponent was unknown since one could not estimate the 
warm H2 column densities. The Infrared Space Observa- 



* ISO is an ESA project with instruments funded by ESA 
Member States (especially the PI countries: France, Germany, 
the Netherlands and the United Kingdom) and with the par- 
ticipation of ISAS and NASA. 



tory (ISO; Kessler et al. 1996), has allowed us, for the 
hrst time, to measure directly the total column density of 
warm gas by observing pure-rotational lines of H2. These 
trace gas with temperatures of a few hundreds Kelvin. ISO 
has also allowed us to study the H 2 ortho-to-para ratio 
(OTPR), which can help determine the possible heating 



2 



N.J. Rodriguez- Fernandez et al.: Non-LTE OTPR in two GC clouds 



mechanism and the origin of this molecule. Before ISO, the 
H2 OTPR had been studied in regions with temperatures 
of ~ 2000 K, using the vibrational lines. In such shock- 
excited sources, one measures an OTPR of ~ 3 (Smith et 
al. 1997), which is the local thermodynamical equilibrium 



(LTE) value for T> 200 K. In contrast, for regions heated 
mainly by ultraviolet (UV) radiation (Photodissociation 
regions [PDRs]), the vibrational lines give OTPRs in the 
range of 1.2-2 (see e.g. Chrysostomou et al. 1993). How- 



ever, these low OTPRs might not be a consequence of an 
actual non-LTE ortho-to-para abundances ratio but a re- 
sult of optical depth effects in the flu oresc ence-pumping 
of the ortho-H 2 (Sternberg & Neufeld |1999| ). Using these 
considerations, one can explain why the PDR in S140 ex- 
hibits an OTPR ~ 2 in the vibrational states but 3 in the 
lowest rotational levels. 

There are two cases of non-equilibrium OTPR mea- 
sured from the pure-rotational lines: the shock excited 



source HH54 (Neufeld et al. |1998[ ) and the PDR associ- 
ated with the reflection nebula NGC 7023 (Fuente et al. 
1999). The first case has been explained using the shocks 
model of Timmermann (1998), which involves transient 



heating by low velocity shocks. To explain the non-LTE 
OTPR in NGC 7023, it was necessary to invoke a dynamic 
dissociation front. 

To investigate the thermal balance of the GC clouds we 
have s electe d 18 clouds from the sampl es of Huttemeister 
et al. ( 1993 ) and Martin-Pintado et al. ( 1997 ) and we have 
observed them with the ISO satellite. In this paper we 
present H 2 observations toward two sources which show 
similar characteristi cs (al so shown in the NH3 studies of 
Huttemeister et al. 1993 ), indicating that their heating 
mechanisms are also very similar. In particular, they show 
a non-LTE OTPR. The detection of OTPRs out of equi- 
librium in the GC clouds gives us new insights into the 
heating mechanism, since the gas must be heated to sev- 
eral hundreds K almost without changing the OTPR of 
cold gas. 

In Sects. 2 and 3, we present observations and results, 
respectively, and in Sect. 4 we discuss the possible heating 
mechanism and the origin of the non-equilibrium OTPR. 



2. Observations and data reduction 

2.1. ISO observations 

We observed the pure-rotational H 2 lines S(0), S(l), S(2), 
and S(3) with the Short Wavelength Spectrometer (SWS; 
de Graauw et al. 1996) on board ISO toward two molec- 
ular clouds. The galactic coordinates and ISO beam sizes 
are given in Table 1. These sources are among the far- 
thest from the dynamical center of the galaxy in our sam- 
ple. M+3.06+0.34 [a(2000)= 17 h 51 m 26!4, 5(2000) = - 
26°08'29'.'4' 



(Bally et al. 1988). The observations were made during 
orbits 313 (S(0) and S(3) lines), and 467 (S(l) and S(2) 
lines). The wavelength bands were scanned in the SWS02 
mode with a typical on-target time of 100 s. The spectral 
resolution (A/ AA) of this mode is ~ 1000-2000 correspond- 
ing to a velocity resolution of ~ 150-300 km s _1 . All the 
lines have broader profiles than those expected for a point 
source by a factor 1.3-1.5, indicating that the sources are 
extended in the direction perpe ndicular to the slit (see 
Valentijn & Van der Werf 1999 ). The flux calibration is 
believed to be accurate to 30%, 20%, 25%, and 25% for 
the S(0), S(l), S(2), and S(3) lines, respectively (Salama 
et al. 1997). Data reduction was carried out with version 
6 of the SWS Interactive Analysis at the ISO Spectrome- 
ter Data Center at MPE. Further analysis has been made 
using the ISAP [j] software package. All lines have been 
rebinned to one fifth of the spectral resolution of the in- 
strument. Fig. [l^i-b shows the spectra, and the observed 
parameters are given in Table 1. The errors in the radial 
velocities of the H 2 lines listed in this table have been esti- 
mated from the Gaussian fits. The wavelength calibration 
uncertainties, expressed in velocities, are typically of 20-40 
km s~ x for A > 12/im and ~ 25-60 km s _1 for A < 12/im 
(Valentijn et al. 1996). Thus, the calibration uncertainties 
usually dominate the global error in the radial velocities. 
When one takes into account the errors from the Gaus- 
sian fits and the wavelength calibration uncertainties, the 
central velocities of the H 2 lines are in agreement with 
those measured from the CO lines (section 2.2). It is note- 
worthy that, the higher the signal-to-noise ratio of the H 2 
lines (S(l) lines), the better the agreement of the H 2 radial 
velocities with those of CO. 

We also present Long Wavelength Spectrometer (LWS; 
Clegg et al. 1996; Swinyard et al. 1996] ) observations of 
these sources in grating mode (43-196.7 /mi, A/AA ~ 200). 
Fig. |l| e-f shows the LWS spectra. The spectral resolution 
was 0.29 /mi for the 43-93 /mi range and 0.6 /mi for the 80- 
196 /mi range. The LWS aperture was ~ 80" x 80" . The 
roll angle, which gives the orientation of the apertures, 
was 90°± 2°for both the SWS and the LWS observations. 
Data were taken during orbits 315 and 318 and processed 
through the LWS Pipeline Version 7. The individual de- 
tector scans were calibrated to within 10% of each other, 
based on overlapping detectors. Post-pipeline analysis (in- 
cluding shifting the different detectors using dark currents 
and dcfringing) was performed with ISAP. 

2.2. IRAM 30-m observations 

The J=l-0 line of 13 CO and C 18 and the J=2-l line of 
C 18 were observed simultaneously with the IRAM 30-m 



is located in the "Clump 2" complex (Stark & 
Bania |l986|) , while M+l. 56-0.30 [a(2000)= 17 h 50 m 26!6, 
(5(2000) = -27°45'29'.'5] belongs to the "Z=l?5-complex" 



The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and Data 
Centers. Contributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 
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telescope (Pico Veleta, Spain) in May 1997. We used two 
SIS receivers at 3 and 1.3 mm connected to two 512 x 1 
MHz channel filter banks. This configuration provided a 
velocity resolution of 2.7 and 1.4 km s~ x for the J=l- 
and J=2-l lines respectively. Typical system tempera- 
tures were ~ 250 K for the J=l-0 line and ~ 500 K for 
the J=2-l line. The receivers were tuned to single side 
band with rejections always larger than 10 dB that were 
checked against standard calibration sources. The beam 
size of the 30-m telescope was 22" and 11 at 3 and 1.3 
mm respectively. Pointing and focus were monitored regu- 
larly. Pointing corrections were always found to be smaller 
than 3 . Calibration of the data was made by observing a 
hot and cold loads with known temperatures, and the line 
intensities were converted to main beam brightness tem- 
perature, Tmb; using main beam efficiencies of 0.74 and 
0.48 at 3 and 1.3 mm respectively. The spectra are shown 
in Fig [j]c-d and the observed parameters as derived from 
Gaussian fits are listed in Table 2. 



3. Analysis 

In Fig. U we show the H2 rotational diagrams for the two 
sources. The open squares correspond to the column densi- 
ties as measured with ISO, without any correction for the 
different apertures in the different lines and for the dust 
extinction. The rotational diagrams for the two sources 
show a zig-zag distribution since the column densities in 
the ortho-H2 levels J=3 and 5 are lower than those ex- 
pected from the para-Hi2 levels for the LTE OTPR. For the 
typical temperatures involved in these transitions (~ 200 
K), the LTE OTPR is ~ 3. Similar rotational diagrams 
derived from the H2 pure-rotational lines h ave b een pre- 
viously found in HH54 by Neu feld et al. (199S) and in 
NGC7023 by Fuente et al. ( [19991) . For these sources where 
extinction is known to be low, the immediate conclusion 
was that the OTPR was not in LTE. 

The H2 emission has been detected in all sources of our 
sample indicating that the H2 emission in the GC must 
be rel atively widespread and extended (Martin-Pintado 
1999b). This is also suggested from the measured 
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linewidths of the H2 lines (see Section 2). Anyhow, even in 
the extreme case that the H2 emission were point-like, the 
corrections for the different apertures would be small and 
would not affect substantially the conclusions about the 
OTPR. For a point-like source the S(0) line will be more 
diluted than the S(l) and S(2) lines because of the larger 
beam (20" x 27" instead of 14" x 27"). The situation for 
the S(3) line will be the opposite since the aperture at this 
wavelength is 14 x 20" . Therefore, in this limit case, the 
column densities in the level 3=2 (derived from the S(0) 
line) averaged in a beam of 14" x 27 would be larger by 
a factor of 1.4, while on the opposite, the beam-averaged 
column density in the J=5 level would be smaller by a 
factor of 1.4. Hence, the correction for different apertures, 



cannot explain the zig-zag distribution in the rotational 
diagram. 

A more critical correction is that for the extinction 
produced by the foreground material. As described by 
Martin-Pintado et al. ( |1999b[ ) the weakness of the S(3) 
line in the GC clouds should be due to the extinction pro- 
duced by the silicate feature at 9.7 /im in the foreground 
dust clouds. In clouds with a LTE H 2 OTPR, one can use 
the intensity of the S(3) line to estimate the visual ex- 
tinction once the relative value for the opacity at 9.7 /im 
to the 0.55 /xm opacity is known. One could, in principle, 
apply corrections for increasing extinctions until the col- 
umn density in the J=5 level is consistent with the column 
densities derived for other levels, i.e., until the rotational 
plot is a straight line (in the case of a Boltzmann distribu- 
tion with one source temperature) or a smooth curve (in 
the case of a temperature gradient). In clouds with a non- 
equilibrium OTPR one could use an equivalent method 
using only ortho-H2 levels, but obviously more than two 
levels are needed. The effect of foreground extinction on 
the rotational diagram is illustrated in Fig. 0, where the 
observed fluxes have been corrected for 30 (filled triangles) 
and 60 mag (filled circles) of visual extinction, using the 
extinction law of Draine & Lee (1984). Visual extinctions 



larger than 60 mag are needed for consistency between 
the S(l) and S(3) line intensities and a LTE OTPR. In 
this case, the curvature of the rotational plots suggests 
the presence of a large temperature gradient in the H2 
emitting region. To constrain the visual extinction toward 
these sources, in the following sections we will estimate 
the total column densities of dust and gas from measure- 
ments of the continuum dust emission, 13 CO, and C 18 
with a similar resolution to that of the SWS aperture. 



3.1. H2 column densities from C 18 O and 
observations 



13 



CO 



We applied the Large Velocity Gradient (LVG) approxi- 
mation to our data, to derive the physical conditions and 
the column densities of molecular gas from the emission 
of the J=2-l line of C 18 and the J=l-0 lines of C 18 
and 13 CO. The lines toward the two sources show com- 
plex profiles with two velocity components. From the line 
intensity ratios one can see that these components have 
slightly different physical conditions. The C 18 J=2-l to 
J=l-0 line ratio is 1.0-1.4 in M+3.06+0.34 and cannot 
be determined for the other source. The J=l-0 13 CO to 
C 18 ratio ranges between 5 and 14 in M+3.06+0.34 and 
is > 7 in M+l.56-0.30. To within a factor of 2, these val- 
ues are in agreement with the typical isotopic abundances 
found in the GC for carbon and oxygen (see Wilson & 
Matteucci 1994) indicating that the 13 CO lines are op- 
tically thin. From the C 18 J=2-l to J=l-0 ratio we de- 
rive for M+3.06+0.34 the H2 densities given in Table 3 
for two cases: high kinetic temperature (Xk=100 K) and 
low kinetic temperature (Tk=20 K). For those H2 densi- 
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ties we have constrained the total column densities using 
the 13 CO line intensities. When the C 18 lines were not 
detected the range of possible 13 CO(1-0) column densi- 
ties was obtained by changing the H2 density between 10 3 
and 10 4 cm" 3 for T K =100 K and between 10 3 5 and 10 45 



cm" 3 for Tk= 20 K (see Hiittemeister et al. 1998). In the 



case of cold gas and even higher H 2 densities, the 13 CO 
column densities will increase only in a factor of 1.3 since 
for low temperatures and densities > 10 4 the J=l-0 tran- 
sition of 13 CO is thermalized. The H2 column densities, 
./Vh 2 , in Table 3 have been derived from the 13 CO column 
density and a fractional abundance with respect to H2 of 
5 10 -6 . They are typically of a few 10 22 cm~ 2 , in good 
agreement with the values given by Hiittemeister et al. 
( 1998 ). With these column densities, we have derived the 
total visual extinction, A v , using the standard conversion 
factor: A r n 2 ( cm~ 2 )=j4 v (mag) x 10 21 . Thus the extinc- 
tions toward the two GC sources studied in this paper are 
typically of 15-20 magnitudes. 



3.2. Dust column densities and temperatures 

From the LWS data we can make a direct estimate of the 
dust temperature and the dust column densities toward 
both sources. Though the aperture of the LWS is larger 
than that of the SWS, the dust emission in the GC is 



relatively smooth (Odenwald & Fazio 1984 ) and one does 
not expect large variations within the LWS aperture. The 
spectra for the two sources have very similar shapes with 
the maximum of the emission at ~ 100 ^m, indicating 
that the bulk of the dust is relatively cold with temper- 
atures below 30 K, in agreement with previous estimates 
(Odenwald & Fazio |1984| Gautier et al. |1984|) . 

The data cannot be fitted with only one gray body. 
For simplicity, we have considered a model with two gray 
bodies of temperatures T\ and The total flux, S\, is 
given by: 



Sa = 0[B(T 1 ,A)(1- 



-(l-/)r(Ah 



-B(T 2 ,\)(l- 



-frW 



)](1) 



where fi is the solid angle of the continuum source, B(T) 
is the Planck function, / is the fraction of the opacity 
due to the warmer component (T2), and r(A) is the total 
opacity at wavelength A. In this model, the ratio of the 
visual extinction, A v , to the total optical depth at 30 fim 
is taken from the Draine & Lee (1984) extinction law and 
the opacity for A > 30 is given by: 



r(A) = 0.014A v (30/im/A) a 



(2) 



where a is the spectral index of the dust emission. In ac- 
cordance with previous estimates for the envelope of Sgr 



B2 (Martin-Pintado et al. |1990J) and for the GC back- 
ground of the cold core GCM 0.25+0.11 (Lis & Menten 
1998| ), we have taken a ~ 1. We have assumed extended 
emission (£1 = SIlws) and then we have fitted the con- 
tinuum spectrum with /, A v , Ti, and Ti as free parame- 
ters. As an example, we show in Fig. || the best fit to the 



LWS spectra towards M+l.56-0.30 obtained with A v =40, 
Ti=15 K, T 2 =27 K and /=0.1. 

Table 4 lists the results of the parameters for the best 
fits for the two sources. The visual extinctions derived for 
the two sources are 30 and 40 mag. These values are in 
agreement, to within a factor of 2, with those derived from 
the CO data. 

The dust emission is dominated by the cool (T ~ 15 
K) component (r Vl ~ (1 — /)t v ), while the slightly warmer 
component (T ~ 30 K) contributes only 10%-20% to the 
total optical depth (r V2 ~ /r v ). We can also fit the spec- 
tra with larger spectral indexes by increasing the dust col- 
umn densities. For instance, an spectral index of 1.5 will 
increase the visual extinction to 50-100 mag. These high 
values of A v are very unlikely since they are almost one 
order of magnitude higher than the estimates made from 
CO (see Table 3). 

Since the extinction derived from CO and the contin- 
uum accounts for the total gas and dust along the line of 
sight, they must represent an upper limit to the extinction 
to the H2 emitting region. Considering the uncertainties 
introduced by the unknown spectral index and the many 
free parameters in the dust column density determina- 
tion, in the following discussion we will assume that upper 
limits to the visual extinction of the H2 emitting region 
are those derived from the CO emission, namely, 16 mag- 
nitudes for M+3.06+0.34 and 20 mag for M+l.56-0.30. 
These values are within a factor of two of estimates ob- 
tained from the total dust column density. 



3.3. Warm Hi: ortho-to-para ratio and column densities. 

As discussed at the beginning of Sect. 3 the H2 OTPR 
depends on the correction for extinction. In the previ- 
ous sections we have estimated the extinction for the two 
clouds and Fig |J shows the H 2 rotational diagrams for 
M+3.06+0.34 and M+l.56-0.30 corrected for the esti- 
mated extinctions. The error bars take into account the 
errors in the Gaussian fits of the lines and the calibration 
uncertainties. From these data, we derive an ortho rota- 
tional temperature, T Q , from the ortho-H2 levels J=3 and 
J=5. In the same way, one can define a para rotational 
temperature, T p , derived from the para-H2 levels J=2 and 
J=4, and an ortho-para temperature, T op , derived from the 
ortho level J=3 and the para level J=2. These tempera- 
tures are listed in Table 5. As we see, T p is ~ 250 K for 
both sources while T D is slightly higher (~ 270 K) indi- 
cating the presence of a moderate temperature gradient. 
This effect is more definite in other sources of our sample, 
where the S(4) and S(5) lines, which trace clearly higher 
te mperat ures, have also been observed (Martin-Pintado et 
al. 1999h| ). For the present sample, T op is ~ 160 K, much 
smaller than T p and T Q indicating a non-LTE OTPR. In 



5 



terms of these temperatures, the OTPR measured from 
our data will be given by: 



OTPR = OTPR LTE (Tp) exp(i- - -L) 

-*■ p op 



(3) 



where OTPR LTE (T) is the LTE OTPR at temperature T. 
As mentioned before, OTPR LTE is - 3 for T > 200 K. 
Using Eq. (||) , one finds an OTPR of ~ 1 for both sources 
(see Table 5). Increasing the extinction will make the H2 
OTPR closer to the equilibrium value, however extinc- 
tions > 70 mag will be required to give an LTE OTPR 
ratio. Such large visual extinctions are very unlikely from 
the molecular line and continuum data discussed in the 
previous sections. We therefore conclude that for the two 
sources the H2 OTPR is not in equilibrium. Since the es- 
timated error is ~ 0.4, we can take ~ 1.4 as a conservative 
upper limit for the OTPR in these two sources. 

Extrapolating the populations in the 3—2 and J=3 lev- 
els to the J=0 and J=l levels , respectively, as two different 
species at temperature T p , one finds that the total column 
densities of warm H2 are ~ 2 10 21 cm~ 2 . This must be 
considered as a lower limit to the actual warm H2 column 
density since the populations of the lowest levels ( J=0 and 
J=l) can be increased significantly by colder, though still 
warm (~ 100 K) gas. Of course, if extinction is higher 
column densities will also increase. This implies that the 
measured ratio of warm H2 to cold gas traced by CO is at 
least 15%. 

High gas kinetic temperatures in these two clouds 
are known to be present from the NH3 observations of 
Huttemeister et al. ( [1993 ). The rotational temperatures 
derived from the (4,4) and the (5,5) metastable inversion 
lines of NH3 are in good agreement with the temperatures 
derived in this paper using the lowest H2 pure-rotational 
lines. Extrapolating the populations in the (4,4) and the 
(5,5) NH3 levels to lower levels with the rotational tem- 
perature derived for each source by Huttemeister et al. 
one finds a column density of warm NH3 of ~ 7 10 



cm -2 in both sources. Taking into account the warm H2 
column densities given above, we find a NH3 abundance 
of (2-4) 10 -7 , similar to the value obtained by Martin- 
Pintado et al. ( |1999a[ ) in the expanding shells of the enve- 
lope of Sgr B2. A similar abundance is also obtained when 
we compare the colu mn d ensities of cold (~ 20 K) NH3 
(Huttemeister et al. 1993 ) and the H2 column densities 



derived by our 13 CO and C 18 data. 
3.4- Warm dust column densities 

If the gas and dust are coupled, one expects that the dust 
associated with the warm H2 component would be an in- 
tense continuum emitter in the mid- and far-IR. There is 
no hint of such dust component in our data, as shown in 
Fig. [| where we represent (as a dotted line) the emis- 
sion of a gray body with a temperature of 250 K and the 
size of the SWS aperture attenuated by the total column 



density of the cold component. The equivalent H2 column 
density of warm dust used to simulate the emission in Fig. 
H is only 5 10 18 cm~ 2 . Even this small column density 
should have been detected. Hence, we can rule out a dust 
component coupled to the warm gas with a column den- 
sity larger than 2 10 -3 times that of the warm H2. On the 
other hand, the comparison of CO emission with the cold 
dust emission shows agreement with the standard gas-to- 
dust ratio within a factor of two. 



4. Discussion 

4-1- Heating of the warm component 

The large column densities of warm H2 and the low col- 
umn densities of associated warm dust require a heating 
mechanism that heats selectively the gas maintaining the 
dust at much lower temperatures. A PDR with an inci- 
dent far-ultraviolet (FUV) flux Go of ~ 100 (measured in 
units of 1.6 10~ 3 ergs cm -2 s _1 ) can heat the gas via 
photoelectric effect in the grains to temperatures of 100- 
200 K in the external layers of the cloud without heating 
the dust to temperatures above 30 K (see Hollenbach et 
al. 1991). However, the large gas phase NH3 abundance, 
as derived in Sect. 3.3, is not possible in such a PDR sce- 
nario. The evaporation temperature of NH3 is ~ 75 K, 



therefore it cannot be evaporated from grain mantles at 
only 30 K. Even in the case that evaporation occurs, the 
UV radiation that heats the dust would destroy the fragile 
NH 3 molecule. This is the behavior found in NGC 7023 
where the NH3 abundance is ~ 10~ 8 in the well shielded 
region and decreases by more than a factor of 30 towards 
the region where the UV radiation increases and the dust 
temperature is ~ 70 K (Fuente et al. 1990). 



Shocks have been invoked as an imp ortant heating 
mechanism for the G C clouds (Wilson et al. 1982 Martin- 
Pintado et al. 



1997 



1998|) . In fact, 



Huttemeister et al. 
M+l.56-0.30 belon gs to the u l =l?5-complex" , where 
Huttemeister et al. ( 1998 ) derived the highest SiO abun- 
dance within their sample, while the CS abundance (which 
traces all dense gas, not just the part that has been sub- 
jected to shocks) is not enhanced. They interpreted the 
SiO enhancement to be produced by large scale dynamic 
effects, proposing that in this complex, gas sprayed from 
the intersection of the x% and X2 orbits is crashing into ma- 
terial that is still on X\ orbits in the context of a bar mor- 
phology. In our sample, M+1.56— 0.30 is also the source 
with the highest SiO to CS ratio. Furthermore, Dahmen 



et al. (1997), studying the HNCO emission in this region, 
found evidence for collisional excitation by shocks. 

On the other hand, M+3.06+0.34 is located close to 
one of the CS cores detected by Stark & Bania ( |1986[ ) in 
the "Clump 2" complex. These dense cores are gravita- 
tionally bound but most of the CO is emitted from the 
l ower density gas, not bound to the cores. Stark & Bania 
( 1986| ) suggested that this material is the result of tidal 
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stripping of the cores. It is definite that shocks can play a 
role to explain the large column densities of warm H2 and 
the relatively large abundances of NH3 , as well as the high 
kinetic temperatures in these clouds. In addition, transient 
heating by shock waves provides a natural explanation for 
H2 OTPRs out of equilibrium. 

We have compared the results of the model calcula- 



10 5 yr. It is worth-noting that recent models for ortho- 



tions for slow shocks by Timmermann (199S) with our H2 



data. Interpolating the H2 line strengths predicted by the 
model for a preshock OTPR=l as a function of density, 
we found that a shock with velocity of 10 km s _1 and 
a preshock H2 density of ~ 2 10 5 cm -3 reproduces the 
observed line intensities. The results are displayed in Fig. 
|^ in the form of rotational plots. Open squares are the 
predicted column densities, while filled circles are the val- 
ues derived from observations after correcting for extinc- 
tion. Though the observed flux in the S(0) seems to be 
slightly larger than in the model, the agreement is excel- 
lent, and calibration errors can account for the discrep- 
ancies. The preshock density seems somewhat high but it 
is plausible since the S(3) line is apparently thermalized, 
which implies a lower limit to the H2 density of ~ 10 4 
cm -3 . In any event, the widespread distribution of the 



HCN emission (Jackson et al. 1996) shows that densities 



of ~ 10 5 cm 3 are common in the GC. 



4-2. The ortho-to-para ratio 

The main processes that affect the OTPR of H2 are proton 
exchange collisions with H + and reactive collisions with H 
atoms. Ortho-para conversion in grain surfaces is thought 
to be less efficient. The rate coefficient for the proton ex- 
change reaction 



H 2 (ortho) + H+ ^ H 2 (para) + H+ + 170.5 K 
3 10- 10 



(4) 



is 



cm 3 s (Gerlich 1990). The analogous re- 
+ 



actions with Hj and HsO + may also occur at a similar 



rate (see e.g. Le Bourlot et al. 1999| ). This rate gives an 



ortho-para conversion timescale, r conv , of ~ 100/n( + j yr, 
where n( + ) represent the density of H + , H^ or H30 + in 
cm~ 3 . One should note that the actual conversion time 
can be a factor of 10 larger than T conv (see Flower and 
Watt 1984). 

The rate coefficient for the reactive collisions with H 
atoms 



H 2 (ortho) + H # H 2 (para) + H + 170.5 K 



(5) 



is 



10 -n e -(3900/T) cm 3 s" 1 (see e.g. Le Bourlot et al. 
1999] ). Due to the high activation barrier of this reaction 
(3900 K) , in cold and dense molecular clouds the dominant 
process will be proton exchange collisions. This is also true 
for low velocity shocks of ~ 10 km s _1 since the maximum 
temperature achieved in the post-shock region is only ~ 
300 K. If the H+ and H+ densities (n(H+), n(H+)) in the 
postshock region of a 10 km s _1 shock were as high as 
~ 10~ 3 cm~ 3 (see Timmermann [l998 ) , r conv would be 



para conversion in shocks by Wilgenbus et al. (2000) find 



much lower H + and H3 densities in the postshock region. 
In this case, the timescale for ortho-para conversion would 
be > 10 5 yr. On the other hand, the time needed for the 
passage of the proposed 10 km s _1 velocity shock, from 
the point where the neutral gas starts to heat up to the 
point where the gas has reached interstellar temperatures 
again, is <10 4 yr (see Timmermann 1998). However the 



timescales in which the neutral gas is at high temperatures 
are much shorter. Hence, if the initial OTPR was lower 
than 3, the heating-cooling of the gas is too fast for the 
OTPR to reach the equilibrium at the temperatures of the 
shocked material. 

Shocks with velocities > 20 km s _1 heat the gas to 
temperatures > 700 K. Then, collisions with H would 
be the main ortho-to-para conversion mechanism, and 
indeed, the ortho-para conversion timescale would be 
low enough to obtain at least some conversion in the 
shock timescale as in the source HH54 (Neufeld et al. 
1998). However, the lines ratios in M+1.56— 0.30 and 



M+3.06+0.34 cannot be explained with a preshock OTPR 
of < 1 and a shock with velocity >10 km s _1 . There- 
fore, the observed OTPR in these clouds must be approx- 
imately the preshock OTPR. This conclusion is indepen- 
dent of any shock model since the low temperatures in- 
volved by a 10 km s^ 1 shock are not sufficient for the 
H2-H reactive collisions to be effective and, even for the 
largest predicted H + and H^ abundances, the proton ex- 
change reactions are not fast enough to give ortho-para 
conversion in the shock timescale. 

If the OTPR of the preshock gas was in equilibrium at 
the gas temperature, the temperature should be ~ 80 K. 
In this case, the preshock gas should have been already 
heated before the shock front compresses and heats the 
gas to 250 K. However, there is no strong reason to be- 
lieve that the preshock OTPR should be in equilibrium 
at the preshock temperature. The H2 molecule is formed 
mainly on the grain surfaces by a highly exothermic reac- 
tion. Thus, if it is rapidly ejected to gas phase the OTPR 
will be the typical OTPR at high temperature, i.e., 3. On 
the other hand, if it is not evaporated immediately from 
the grain there will be ortho-to-para conversion by colli- 
sions with radicals, impurities or defects and the OTPR 
could reach the equilibrium value at ~ 30 K (dust tem- 
perature) of ~0.01. In our case, the preshock OTPR of 
~ 1 suggests that the H 2 molecules were ejected from the 
grains with OTPR >1. Afterwards, this ratio could de- 
crease due to proton exchange processes. 

The equilibrium proton abundance in dense (n(H2) 
~ 10 5 cm~ 3 ) clouds, where photoprocesses are not im- 
portant, depends mainly on the ionization by cosmic rays 
and on charge exchange reactions with neutral molecules. 
Modeling the c hemistry of dense PDRs, Sternberg and 
Dalgarno (|199ED found n(H + ) of ~ 10 5 cm 3 in the well 
UV-shielded region for a cosmic ray ionization rate (£) of 
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5 10 — s _1 , implying r conv ~ 10 7 yr. A similar timescale 
is obtained for proton exchange collisions with H^~. The 
density of H 3 + could reach 10~ 4 cm~ 3 and thus T conv 
could decrease by a factor of 10. Nevertheless, the actual 
time to reach the LTE OTPR would be longer. Flower & 
Watt (1984) have studied the temporal evolution of the 
OTPR in molecular clouds. Using the same rate coeffi- 
cient as above for the proton exchange process, they have 
shown that for H + densities ^] of 10~ 4 — 1CP 5 cm -3 the 
actual time needed for an OTPR=3 to be in equilibrium 
at 30 K (similar to the observed cold component in the 
GC clouds) is ~ 10 7 — 10 s yr. In particular, if n(H + ) (or 
n(H 3 0+)) is - 10~ 4 cm~ 3 , then - 5 10 6 yr will be 
needed to have an OTPR=l. Assuming that the H2 was 
ejected to gas phase after formation with an OTPR <3, 
the shock front reached the cloud approximately 10 6 yr 
after the formation of the H2 molecules, since this is the 
time needed for an OTPR^3 to descend to ~1 in a dense 
molecular cloud. 



5. Conclusions 

We have presented ISO SWS observations of the S(0), 
S(l), S(2), and S(3) pure-rotational lines of H 2 and LWS 
observations of the dust continuum and IRAM-30m 13 CO 
and C 18 observations toward the GC molecular clouds 
M+l. 56-0.30 and M+3. 06+0.34 . Using the CO data and 
dust column densities from the LWS spectra we estimate 
~ 20 mag of visual extinction toward these sources. The 
two estimates of the extinction agree within a factor of 2 
for the standard gas-to-dust conversion factor. According 
to the two c omponents scenario proposed by Huttemeister 
et al. ( 1993 ), the low- J CO emission arises from the cold 
gas component and is coupled to the dust at a tempera- 
ture of < 30 K. The warm component (T ~ 250 K) column 
density observed directly in H2 is, at least, ~ 15 % of the 
cold one and would have very little warm dust associated 
with it. From the LWS spectra we set a conservative upper 
limit to the warm gas column density associated with the 
warm dust of 2 10~ 3 times that of the warm H2 column 
density. 

After correcting for the dust extinction, we derive an 
OTPR of 1.0 ± 0.4, which is far from the LTE value ex- 
pected for the gas temperatures of 250 K. We have also 
compared the warm H2 column densities to the NH3 ob- 
servations by Huttemeister et al. (1993), and derived NH3 
abundances of ~ 2 10 -7 similar to those in the cold com- 



ponent (Huttemeister et al. 1993 ). 

The low dust temperatures, the high NH3 abundances, 
the large CO linewidths, the non-LTE H 2 OTPR, in addi- 
tion to the high gas temperatures suggest that the warm 
gas component is heated by low velocity shocks with 
speeds of ~ 10 km s _1 . To explain the OTPR we pro- 



2 In their model these proton densities were obtained with 
£ = 10 -17 — 10 -18 s -:L using a simplified chemical network. 



pose the following scenario. H 2 is formed in the grain sur- 
faces and ejected to gas phase with OTPR<3. After ~ 10 6 
years, the time required to reach the preshock OTPR=l, 
a low velocity shock heated the gas to the observed tem- 
peratures of 250 K, but the OTPR was almost unaltered 
because the timescale for the passage of such a shock is 
much shorter than the ortho-to-para conversion timescale. 
Taking into account the shock timescale this occurred less 
than 10 4 yr ago. It is interesting to note that the timescale 
of the cloud's galactic rotation period is also ~ 10 6 years. 
This fact suggests that the origin of the shocks can be 
related to large scale dynamics of the GC region. 
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Table 3. Derived physical conditions for the 13 CO and 
C 18 lines for T K =20 K and T K =100 K 



Source 


T K 


log(n ff2 ) 


N13 co " 


Nh 2 




K 


log( cm" 3 ) 


10 16 cm" 2 


10 21 cm" 2 


M+3.06+0.34 


100 


3-3.8 


5.8-8 


12-16 




20 


>3.3 


<8 


<16 


M+l.56-0.30 


100 


3-4" 


2.8-10.2 


5.7-20 




20 


3.5-4.5° 


2.8-8.1 


5.6-16.2 



used, 

-4 



following 
3 when 



When C O is not detected, we have 
Hiittemeister et al. (1998), rtu 2 ~ 10 3 ' cm 
T K =100 K and, n H2 ~ IO 3 ' 5 - 4 ' 5 cm -3 when T K =20 K. 
Contains the contribution from both velocity components. 
Dispersion in the column densities are due to errors from 
the Gaussians fits to the spectra. 



Table 4. Parameters of the best fits to the LWS spectra 
with two gray bodies assuming a=l and Q = fluws- Tem- 
peratures, ratio of the opacity in the warmer component 
to the total opacity, and total visual extinction. Numbers 
in parentheses are la errors of the last significant digit. 



Source 


Ti 


T 2 


/ 


A v 


M+3.06+0.34 


14(4) 


24(2) 


0.2(2) 


30(20) 


M+l.56-0.30 


15(4) 


27(3) 


0.1(1) 


40(20) 
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M+3.06+0.34 




i — i — i — | — r~ 
a) 



M + 1.56-0. 30 
~i — i — I — i — i — i — i — I — i — i — i — r 



b) 



/ v tor (km s ) \ 



-500 x \ 500 
/ v tar (km s" ) \ 
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d) 
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v,, r (km s ) v, sr (km s" ) 
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5 10 ™ - 
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e) 
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Fig. 1. Spectra of the two sources: a-b H2 spectra taken with the SWS. c-d IRAM-30m spectra of the C 18 O(l-0), 
13 CO(1-0), and C 18 0(2-l) lines, e-f LWS full grating spectra. Note the different radial velocities ranges in Fig. la-b 
and Fig. lc-d. 

Table 1. Observational parameters of the H 2 lines. Fluxes in units of 10 -20 W cm~ 2 . Heliocentric velocities (v\ SI ) 
in km s _1 . Numbers in parenthesis are la errors of the last significant digit of the Gaussian fits. 





S(0) 


S(l) 


S(2) 


S(3) 


Wavelength 


28.2^m 


17.0^m 


12.8^m 


9.7^m 


Aperture 


20" x 27" 


14" x 27" 


14" x 27" 


14" x 20" 




Flux v\ BI 


Flux v\ BI 


Flux v\ BI 


Flux v\ BI 


M+3.06+0.34 


6.9(6) -13(7) 


19.6(6) 3(4) 


16(3) -30 (20) 


1.6(5) 80(20) 


M+l.56-0.30 


8.1(9) -80(10) 


19(1) -43(4) 


14(3) -50(20) 


1.1(4) -30(20) 
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Fig. 2. Rotational plots. The results are displayed for three different values of the visual extinction: 0, 30, and 60 mag. 
We can see the typical zigzag distribution of a non-LTE OTPR. Extinctions higher than 60 mag are needed to have 
the smooth characteristic curve of emission arising from gas with an equilibrium OTPR and a temperature gradient. 




Fig. 3. The LWS spectrum of M+1.56— 0.30 . The thick solid line is the best fit with two components of temperatures 
15 K (dashed line) and 27 K (dot-dashed line). The solid line is the total emission for an equivalent column density 
of iV H2 = 5 x 10 18 enr 2 of hot (250 K) dust with Q = 20" x 20" (dotted line) located behind the cold dust. We have 
assumed that the hot component is extincted by the cold component. 
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Fig. 4. Rotational plots after correcting for the most probable extinctions. The slope of the lines is proportional to 
the inverse of the temperature. T p is the rotational temperature between the para-H 2 levels, T between the ortho-H 2 
levels, and T op between the ortho-level J=3 and the para-level J=2. Error bars take into account the calibration 
uncertainties and the errors in the Gaussian fits. 
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Fig. 5. Filled circles and error bars as in Fig. 4 Empty squares are the expected column densities using the model of 
Timmcrmann (1998) with preshock density ~ 2 10 5 cm~ 3 , shock velocity 10 km s _1 , and preshock OTPR=l 
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Table 2. Observational parameters of the C 18 and 13 CO lines. Digits in parentheses are the errors in the last 
significant digits (rms of the Gaussian fits). For C 18 0(2-l) lincwidths were fixed, when detected, to that of the 
C 18 O(l-0) lines. 



Source 


Position 


C ls O(l-0) a 


C ls O(2-l) a 


1,5 CO (1-0) 




C C ° 


Av T MB 


Av Tmb 


Av Tmb 




km s 1 


km s" 1 K 


km s" 1 K 


km s^ 1 K 


M+3.06+0.34 


10.5(4) 


6(1) 0.80(5) 


6(0) 0.8(2) 


10(1) 2.3(2) 




27.9(6) 


8(3) 0.44(5) 


8(0) 0.6(2) 


18(2) 2.6(2) 


M+l.56-0.30 


-51(1) 


<0.14 


<0.20 


19(5) 0.8(4) 




-24.1(6) 


<0.14 


<0.20 


21(1) 1.8(4) 



a Limits are 3cr assuming the same Av than that for the 13 CO. 



Table 5. Para, ortho, and ortho-para rotational temperatures, OTPR and H2 column densities. Numbers in parentheses 
are la errors. 



Source 


T P 


T 


1 op 


OTPR 


Nh 2 




(K) 


(K) 


(K) 




(10 21 cm" 2 ) 


M+3.06+0.34 
M+l.56-0.30 


260(30) 
250(20) 


280(20) 
270(20) 


160(20) 
160(20) 


0.9(0.4) 
1.0(0.4) 


2.6(1.0) 
2.1(0.8) 



